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We calculate fermion-antifermion- "meson" three-point functions in noncompact lattice QED with dynamical 
staggered fermions and use them to extract effective Yukawa couplings. The results are consistent with the 
hypothesis that QED is trivial. 



In strongly coupled lattice QED chiral symme- 
try is spontaneously broken, whereas for small 
bare charges one expects an unbroken chiral sym- 
metry [1]. This observation naturally leads to the 
question if a nontrivial continuum limit can be 
defined at the corresponding critical point. As- 
suming that this is the case one would also like 
to know the nature of the associated continuum 
theory. To investigate these problems several 
groups have performed extensive Monte Carlo si- 
mulations of lattice QED using various techniques 
(see e.g. [2-7]). Yet a consensus has not been re- 
ached (see also [8]). 

From correlation functions of composite 
fermion-antifermion operators one finds at least 
two fermion-antifermion (bound) states ("meson" 
states): a pseudoscalar state (P, the Goldstone 
boson associated with the spontaneous symmetry 
breaking) and a scalar state (S) [5,9,10]. (Effec- 
tive) Yukawa couplings of these states to the fer- 
mions can be extracted from suitable three-point 
functions, whose investigation is the subject of 
this talk. The study should provide additional 
information about the model. In particular, if 
the theory is trivial, the couplings should vanish 
as one approaches the critical point. 

We have performed simulations with dynamical 
staggered fermions using the noncompact formu- 
lation of the gauge field action. The total action 
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is then given by 5 = Sq + Sp with 

1 " ^ _ (2) 

where j3 = 1/e^ (e = bare charge) and 

^^{x) = {-lY^+-+^^-\ (3) 

Ft,r,{x) = Af,(x)+A^(x+fi)-Af,(x+iy)-A^(x).(4:) 

We work on a 12^* lattice choosing periodic boun- 
dary conditions for the gauge field and peri- 
odic (antiperiodic) spatial (temporal) boundary 
conditions for the fermion field XjX- 

Fermionic observables like the fermion propa- 
gator or the three-point functions to be studied 
need gauge fixing. We work with the Landau 
gauge, which can be implemented exactly in 
the noncompact formulation. Since the zero- 
momentum mode of the gauge field does not ave- 
rage to zero in our ensembles, we form blocks 
of, e.g., 20 successive configurations, on which 
this mode is approximately constant, and perform 
our analysis in each block separately taking the 
zero-momentum mode into account explicitly (cf. 
[10]). From these block results we then calculate 
the averages and errors. 
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According to the triviality scenario, the renor- 
malized charge vanishes in the continuum limit. 
Furthermore, the chiral phase transition is well 
described by logarithmically improved mean field 
theory, i.e., the equation of state can be derived 
from an effective linear cr-model [10]. In parti- 
cular, the behaviour of the effective quartic coup- 
ling A is governed by a (one-loop) renormalization 
group equation of the form 



dA 
'dt 



cA^ J 



In a , 



(5) 



where a = (xx) ^nd c > 0. Solving this equation 
In (T I ~ "'^ as a . For a 



one finds that A 



Yukawa coupling g one obtains in leading order 
[11] 



dt 



c' > 0. 



(6) 



Hence one expects on the basis of the triviality 
scenario that g vanishes like 
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ln,7|-i/2. 



(7) 



In order to measure effective Yukawa couplings 
and to compare them with (7) we have calculated 
the three-point function 



G'f\p,q;to,h) = J2' 



2iq-x' — 2ip-x 



X (x(2£, h)Maip- q, 0)x(2f', to))c 



(8) 



Here (• • ■)c denotes the fermion-line connected 
part of the expectation value, the lattice vectors 
x,x' label spatial cubes of size 2^, and Ma(po,t) 
is the standard pseudoscalar (a = P) or scalar 
(a = S) "meson" operator: 



(9) 



with Lpp{y) = i-iy'+y^+y^,cpsiy) = i- in 

our calculations, the spatial fermion momenta 
p,q take the values {p,q) = (0,0), (^63,0), 
(^63,^63) with the lattice size L = 12. The 
spatial "meson" momentum is then given by po = 
p- q. 

In order to define effective Yukawa couplings 
we compare the Monte Carlo data for the three- 
point functions with tree-level formulas derived 



from an effective lattice action describing the in- 
teraction between staggered fermions x,X a 
(pseudo-)scalar field $ via coupling terms of the 
form 



-ffs ^<I'(«)x(«)x(«) (scalar), 
-gpJ2<S>(x)x(x)x(x)(-ir+^^+^' 

X 

(pseudoscalar) . 



(10) 
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In the tree-level formulas for the three-point func- 
tions we replace the free propagators by the full 
propagators calculated in the simulations, in- 
clude the appropriate wave function renormali- 
zations Zp (for the fermion field) and Za (for the 
"mesons") and Fourier transform also the to,ti 
dependences to obtain the following definition 
of the momentum-dependent Yukawa couplings 

9a{p,q)- 

Y.e-^P^'^+^^^'°G^^\p,q;to,h) 
= ga{p,q)ZF{p)~'^ZF{q)~'^Z„{p- q)-2 
(~8) a u\^i(P4-<i4)t 



V3 



X — (M„ (q - p, t)M„ [p - q, 0)), 
V3 



^^"„(cJ)e-'(P-«)^^ (12) 
{^ry4(^)*°e'«^*° 

to 

^ e^'? (-^--^') (X(2£ + cJ, to)x(2f', 0)) } 

X ,x' 

[^r?4(t^)*^e'(f''+'^)*i 
ti 

^ e^^^vi^-^'){^{2x + t,)x{2x' , 0)) }* , 



where Osit) = I, ep{t) = (-1)*, hi^) = m(^), 
Op(uj) = 1. The spatial volume of the lattice is 
denoted by V3 and uj runs over the 8 three- vectors 
with components 0,1. The wave function renor- 
malizations Zp,Za are defined in the usual way 
by comparing the one-particle contributions in 
the full propagators with the corresponding free 
propagators. 
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Figure 1. Effective Yukawa coupiing for tlie Gold- 
stone boson. 

Let us now discuss tlie numerical results. It 
turns out that ga{p,<l) is essentially real: The 
imaginary part is relatively small and noisy. 
Hence we shall ignore it in the following and treat 
(j as a real quantity. In order to obtain a final 
unique answer one would like to send the mo- 
menta to zero. However, due to the antiperi- 
odic temporal boundary conditions for the fer- 
mions we can only give the "mesons" momentum 
zero by choosing p = q. Furthermore, we may 
take p = q = 0. Instead of considering the coup- 
lings at \q4\ = tt/L, the smallest possible value, 
we eventually average (jf„((0, 54), (0, 54)) over q4 
encouraged by the observation that the data for 
this quantity show hardly any dependence on q4. 
For the error we take the largest error of a single 
value. 

This preliminary analysis leads to the data for 
gp plotted in fig. 1 versus | In (rj""'^''^ as sugge- 
sted by eq. (7). They are consistent with a li- 
near extrapolation to zero thereby confirming the 
triviality scenario. The data for gs lead to the 
same conclusion although less convincingly. This 
could be due to the fact that we have neglected 
the fermion-line disconnected contributions to the 
three-point functions, which are very difficult to 
calculate. Fig. 1 also demonstrates that our data 
are still quite far away from the continuum limit: 



The extrapolation has to bridge a considerable 
gap thus leaving room for speculations about non- 
trivial behaviour. 
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